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Q\ , Abstract 

The Earth, Mars, Sun, Jupiter system allows for a sensitive test of the strong equivalence principle 
(SEP) which is qualitatively different from that provided by Lunar Laser Ranging. Using analytic 
and numerical methods we demonstrate that Earth-Mars ranging can provide a useful estimate of 
the SEP parameter 77. Two estimates of the predicted accuracy are derived and quoted, one based on 
conventional covariance analysis, and another (called "modified worst case" analysis) which assumes 
that systematic errors dominate the experiment. If future Mars missions provide ranging measure- 
ments with an accuracy of cr meters, after ten years of ranging the expected accuracy for the SEP 
parameter r\ will be of order (1 — 12) x 10~ 4 er. These ranging measurements will also provide the 
most accurate determination of the mass of Jupiter, independent of the SEP effect test. 

Subject headings: celestial mechanics, stellar dynamics - gravitation - Earth - planets and satel- 
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, lites: Mars - dark matter 
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On 

"q ' 1 Introduction 

cr 

The question, posed long ago by Newton, of the relation between the gravitational and inertial masses of 
the same body continues to be the subject of theoretical and experimental investigations. This question 
arises in most any theory of gravitation. The equality of inertial and passive gravitational masses, often 
^ . stated as the weak equivalence principle (WEP), implies that different neutral massive test bodies will 

have the same acceleration of free fall Aq in an external gravitational field, and therefore in freely falling 
inertial frames the external gravitational field appears only in the form of a tidal interaction (Singe 1960). 
Up to these tidal corrections, freely falling bodies behave as if external gravity were absent (Bertotti & 
Grishchuk 1990). In the construction of the general theory of relativity Einstein went further, postulating 
that not only mechanical laws of motion, but all non-gravitational laws should behave in freely falling 
frames as if gravity were absent. If local gravitational physics is also independent of the more extended 
gravitational environment, we have what is known as the Strong Equivalence Principle (SEP). 

Various experiments have been performed to measure the ratios of gravitational to inertial masses 
of bodies. Experiments on bodies of laboratory dimensions verify the WEP to a fractional precision 
SAq/Aq = 5{m g /m) ~ 10~ n by (Roll, Krotkov & Dicke 1964) and more recently to a precision SAq/Aq ~ 
10" 12 by (Braginsky & Panov 1972; Adelberger et al. 1994). The accuracy of these experiments is 
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sufficiently high to confirm equal strong, weak, and electromagnetic interaction contributions to both 
the passive gravitational and inertial masses of the laboratory bodies. This impressive evidence for 
laboratory size bodies does not, however, carry over to celestial body scales. 

Laboratory size bodies used in the experiments cited above possess a negligible fraction of gravita- 
tional self-energy and therefore such experiments indicate nothing about the equality of gravitational 
self-energy contributions to the inertial and passive gravitational masses of the bodies (Nordtvedt 1968a). 
Interesting results for celestial bodies are obtained if one includes terms of fractional order (f^/mc 2 ), 
where m is the mass of a body B and VLb is its gravitational binding or self-energy: 



(A) (i; 

\ mc J b 2mc z Jv B \x — y\ 



This ratio is typically ~ 10~ 25 for bodies of laboratory sizes, so experimental accuracy of a part in 10 12 
sheds no light on how gravitational self-energy contributes to the inertial and gravitational masses of 
bodies. 

To test the SEP one must utilize planetary-sized extended bodies in which case the ratio (1) is 
considerably higher. Numerically evaluation of the integral of expression (1) for the standard solar 
model (Ulrich 1982) obtains 

-3.52 • 1(T 6 (2a) 



mc 2 J s 

and the analogous value has been obtained for the Earth (Allen 1985): 

n 



mc 2 J E 



4.6 -lO" 10 , (26) 



The development of the parameterized post-Newtonian (PPN) formalism ( Nordtvedt 1968b; Will 
1971; Will & Nordtvedt 1972), allows one to describe within the common framework the motion of 
celestial bodies in external gravitational fields within a wide class of metric theories of gravity. Within 
the accuracy of modern experimental techniques, the PPN formalism becomes useful framework for 
testing the SEP for extended bodies. In that formalism, the ratio of passive gravitational to inertial 
mass is given by (Nordtvedt 1968a,b) 

— y - = 1 + rj 2 ' ( 3 

in which the SEP violation is quantified by the parameter rj. In fully-conservative, Lorentz-invariant 
theories of gravity the SEP parameter is related to the PPN parameters by 

r, = 4/3 - 7 - 3 (4) 



and is more generally related to the complete set of PPN parameters through the relation 

?? = 4/3- 7 -3-y^-ai + i (2a 2 - 2£ - C2). (5) 

A difference between gravitational and inertial masses produces observable perturbations in the 
motion of celestial bodies in the Solar System. By analyzing the effect of a non-zero rj on the dynamics 
of the Earth-Moon system moving in the gravitational field of the Sun, Nordtvedt (1968c) found a 
polarization of the Moon's orbit in the direction of the Sun with amplitude Sr ~ 77C0, where Co is a 
constant of order 13m. We call this effect, generalized to all similar three body situations, the "SEP 
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polarization effect" . The most accurate test of this effect is presently provided by Lunar Laser Ranging 
(LLR) (Williams 1976; Shapiro et al. 1976; Dickey et al. 1989), and in the most recent results (Dickey 
et al. 1994; Williams et al. 1995) the parameter r\ was determined to be 

rj = -0.0005 ± 0.0011 . (6) 

Other tests of SEP violation have been discussed. An experiment employing existing binary pulsar 
data has been proposed by Damour and Schafer (1991). A search for the SEP polarization effect in 
the motion of the Trojan asteroids was suggested in (Nordtvedt 1968a) and carried out by (Orelana & 
Vucetich 1992). Also results are available from numerical experiments with combined processing of LLR, 
spacecraft tracking, planetary radar and Very Long Baseline Interferometer (VLBI) data (Chandler et 
al. 1994). 

It has been observed previously that a measurement of the Sun's gravitational to inertial mass ratio 
can be performed using the Sun-Jupiter-Mars or Sun-Jupiter-Earth system (Nordtvedt 1970; Shapiro 
et al. 1976). This is the first paper from a planed series addressing the above problem. The question 
we would like to answear first is how accurately can we do this ranging experiment? We emphasize 
that the Sun-Mars-Earth-Jupiter system, though governed basically by the same equations of motion 
as Sun-Earth-Moon system, is significantly different physically. For a given value of SEP parameter rj 
the polarization effects on the Earth and Mars orbits are almost two orders of magnitude larger than 
on the lunar orbit. In this work we examine the SEP effect on the Earth-Mars range, which has been 
measured as part of the Mariner 9 and Viking missions. Moreover, future Mars missions, now being 
planned as joint U.S.-Russian endeavours, should yield additional ranging data. 

The dynamics of the four-body Sun-Mars-Earth-Jupiter system in the Solar system barycentric 
inertial frame were considered. The quasi-Newtonian acceleration of the Earth (E) with respect to the 
Sun (S) is straightforwardly calculated to be: 



R SE ym * /l:J k=M,J 
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+ 4(A)c-(A)J E ^§r = A N + A tid + A v (7) 

where /j,* se = /j. s + ^e + V Vs(^) e +Ve(^t) s and fi k = Gm k . The subscripts (M) and (J) indicate 

Mars and Jupiter, respectively. Also Rbc = Rc ~ Rb is the vector from body B to body C and An 
is the Newtonian acceleration term. A^ is the SEP acceleration term, which is of order 1/c 2 . While 
it is not the only term of that order, the other post-Newtonian 1/c 2 terms (suppressed in eq.(7)) do 
not affect the determination of rj until the second post-Newtonian order (~ 1/c 4 ). Finally A tid is the 
Newtonian tidal acceleration term. Note, that A^/An ~ f] ■ 10~ 10 and A tid /A N ~ 7 • 10~ 6 . Given 
that level of accuracy, we ignore the mutual attraction of the two planets, Earth and Mars. The SEP 
acceleration is treated as a perturbation on the restricted three-body problem, and the SEP effect is 
evaluated as an alteration of the planetary Keplerian orbit. 

Using expression (3) and noticing that ixm/Rsm "C [1j/Rsj, we obtain from eq.(7), 

A E -A S & 

* Rse \Rjs Rje] . ( \ Rjs / x 
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Corresponding equations for Mars are obtained by replacing subscript E by M in eqs.(7) and (8). To 
good approximation the SEP acceleration A v has constant magnitude and points in the direction from 
Jupiter to the Sun, and since it depends only on the mass distribution in the Sun, the Earth and Mars 
experience the same perturbing acceleration. The responses of the trajectories of each of these planets 
due to the term A v determines the perturbation in the Earth-Mars range and allows a detection of the 
SEP parameter r/ through a ranging experiment. 

The presence of the acceleration term A v in the equations of motion results in a polarization of the 
orbits of Earth and Mars, exemplifying the planetary SEP effect. We investigate here the accuracy 
with which the parameter r] can be determined through Earth-Mars ranging and approach the problem 
with a series of successive approximations. In Sections II and III the "tidal term" in equation 
(8) is neglected. In Section II the perturbation theory about circular, coplanar reference orbits for 
Earth and Mars is performed. A covariance analysis is carried out to estimate the accuracy to which 
the SEP parameter rj can be determined from a lagre number of Mars ranging measurements, each of 
accuracy a meters. In Section III the calculations of Section II are improved by employing numerical 
integration rather than perturbation theory. The agreement between the two approaches is good, and 
the eccentricity corrections are found to improve the accuracy of the analytic approximation significantly. 
In Section IV the tidal acceleration term A^, is restored, requiring the addition of the mass of Jupiter 
fij to our set of covariance parameters. This mass can be determined more accurately from a few years 
of Earth-Mars ranging than from the Pioneer 10,11 and Voyager 1,2 flybys combined, independent of 
the r\ measurement. In Section V we summarize and suggest further avenues for testing SEP violation. 



2 Perturbation About a Circular Reference Orbit 

Here and in the next section the problem is simplifyed by ignoring the tidal term A^ in equation (8) 
and the corresponding equation for Mars. We examine the effect of the SEP acceleration term A n in (8) 
on the orbits of Earth and Mars by carrying out first-order perturbation theory about the zeroth order 
orbits of Earth and Mars, taken to be circular. Jupiter's orbit is also taken as circular and coplaner 
with Earth and Mars. With these approximations a nine parameter covariance analysis is carried out to 
estimate how precisely Mars ranging can determine the SEP parameter rj. These approximations give a 
standard deviation for r] which closely agrees with the later numerical integration result of Section III. 

A heliocentric reference frame rotating with Jupiter at constant angular frequency uj = \J ps/ R%j is 
assumed. The SEP perturbation is represented by a constant acceleration g„ directed from Jupiter to 
the Sun. Locating Jupiter on the x axis, the following Hamiltonian for both Earth and Mars in polar 
coordinates results: 

H = \{p 2 r + 4) - cope ~ — + 9 v r cosd, (9) 
2 r A r 

where 

Sv^vi—z) Pr = r, Pe = r 2 (9 + u). (10) 

We take the reference orbits for Earth and Mars to be circular Keplerian with r = a and g v = 0. Then 
Pe = \/H-s a an d n = \J ' ns/a %i is the orbital angular frequency. A covariance analysis can be performed 
to estimate the accuracy to which one can determine the SEP parameter through Earth-Mars ranging. 
Earth-Mars range pem is defined as follows: 

Pem = r E + r 2 M - 2r E r M cos(9 M - E ). (11) 
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The variation of bpEM is then: 



fe(t) = E«ft^w, (12) 
fc=l 



where 5pem does not depend on Sco and q is the following vector: 

<! = ( r Eo > r M„ , P£r , PMro > P£0 O j PM6»o > ##0 ~ ^M , PS, 5r?) • (13) 

The partial derivatives are easily computed and for example, 

OpEM 1 



dr E o Pem 



M cos(9e - Om — nst) + -^cos(9e - 9m + n E t) 



a,ECOs(nEt) — 

(14) 



2 v 2 
By definition, the covariance matrix ajk has elements: 

^ 1 dpEM , x dpEM , x n ^ 

i=l a i d 9j d 1k 

where N ranging measurements have been made at times ti,i = 1, N, and have uncertainties Uj. The 
uncertainty in the estimations of the SEP acceleration term g v and parameter rj are then: 



a 9v 



Q \ pj 



mc 2 J s Rj s 



y/i"- 1 )*,*,* ( 16 ) 



with the fractional binding energy (Q/mc 2 ) s given by expression (2a). 

The results obtained are presented in Figure 1, with the thin dashed curve showing the result of 
evaluating equations (15) and (16) for a n assuming N daily range measurements have been taken during 
the mission, each with the same uncertainty a, measured in meters. The initial angles between Earth 
and Jupiter and Mars and Jupiter were taken from the DE242 emphemeris at time 2441272.75, the 
beginning of the Mariner 9 ranging measurements. The quantities he and um are taken to be the mean 
motions of Earth and Mars, cie and their mean distances from the Sun, and the frequency to is taken 
to be the mean motion of Jupiter. The uncertainty in r] first drops very rapidly with time and then 
after a few years approaches the asymptotic behavior ~ iV -1 / 2 . This result gives a lower bound on the 
uncertainty as predicted by conventional covariance analysis. For a mission duration of order ten years, 
the uncertainty behaves as 

a v ~ .0028a/VN (17) 

This result, eq.(17), assumes Gaussian random ranging errors with a white spectral frequency distri- 
bution. But past ranging measurements using the Viking Lander have been dominated by systematic 
error (Chandler et al. 1994). One approach to accounting for systematic error is to multiply the formal 
errors from the covariance matrix by V^V (Nordtvedt 1978). With this approach, the expected error 
decreases rapidly near the beginning of the data interval, but for large N approaches an asymptotic 
value as demonstrated by Fig.l. However, we believe this is overly conservative. A more optimistic 
error estimate would include a realistic description of the time history of the systematic error. For 
example if we knew that the systematic error was a sinusoid of known angular frequency, we could add 
its amplitude and phase as additional variables in the covariance analysis. By contrast the worst-case 
error analysis treats each variable of concern equally by multiplying its computed error by ^/N. 
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A realistic systematic error budget for ranging data to Mars, or for Mercury as considered by a 
group at University of Colorado (Ashby et al. 1995), is not presently available. But it is unlikely that 
we will be so unfortunate that the frequency spectrum of the signal will correspond to the spectrum 
of the systematic error. However, we reduce the upper error bound determined by the s/N multiplier 
(a v = 0.0028<r) by a numerical factor. The ranging experiment proposed by Ashby et al. (1995) for 
Mercury is quite similar to our proposed experiment using Mars. We therefore follow the Colorado 
group and reduce the worst-case error estimate by a factor of three and call the result the modified 
worst-case analysis. This yields an asymptotic value for the error of cr„ = 0.0009<r, in our opinion a 
realistic estimate of the upper bound on the error. 

In the case of the existing Mars ranging derived from the Mariner 9, Viking, and Phobos missions, 
the rms ranging residual referenced to the best-fit Martian orbit was 7.9 m. We computed the covariance 
matrix with assumed daily range measurements for Mariner 9 (actual data interval JD 2441272.750 to 
JD 2441602.504) and Viking (actual data interval JD 2442980.833 to JD 2445286.574). Additionally, 
one ranging measurement from Phobos (actual time JD 2447605.500) was included, although it had 
negligible effect on the result. With a = 7.9 m, a formal error a v = 0.0005 is obtained from the 
covariance matrix. If 7.2 years of Mars ranging is assumed, though not continuous, we reach the 
asymptotic limit of the modified worst-case analysis (as shown by Fig. 1), a realistic error = 0.009 
which is about 17 times the formal error. This is a factor of eight larger than the realistic error set 
by Chandler et al. (1994) from an analysis of the actual combined LLR and Mars ranging data. We 
conclude that the best determination of r\ is provided by the LLR data, but the existing Mars ranging 
can provide an independent solar test with a realistic accuracy interval of 

a v srf 0.0005 - 0.009 (Mariner 9, Viking, Phobos). (18) 

Expression (35) establishes the interval for expected accuracy a v with the lower and upper bounds 
estimated for existed data by conventional covariance analysis and "modified worst case" analysis re- 
spectively. This interval will be narrowed by ongoing upgrades to DSN instrumentation and better 
modelling of the antenna and spacecraft ranging systems (Anderson et al. 1985). Future Mars Orbiter 
and Lander missions are expected to achieve an rms systematic ranging error between 0.5 and 1.0 m. 
Then after a few years of ranging, the realistic error on rj should fall to around 

a v (0.00006 - 0.0011) a (Future Mars missions). (19) 

3 Numerical Integration Without the Tidal Term A t id- 

To obtain a more accurate estimation of and to check the results of Section II, a numerical integration 
of the Sun-Earth-Mars-Jupiter system in heliocentric coordinates is performed. Equation (8) without 
the A t id term was used for Earth and Mars (with the interchanging of the subscripts {E) — ► (M)). 
The equation of motion for Jupiter did not include a SEP term or Newtonian perturbation from other 
planets. The DE242 ephemeris was again used for the initial conditions, and the same 9 parameters of 
equation (30) were used in the covariance analysis. 

Assuming frequent Earth-Mars range measurements as in Section II, the results are shown in Figure 
1 (thick dashed curve). The results agree fairly well with the analytical ones of Section II (thin dashed 
curve). For ten years of observations we obtain 

a v ~ .0027a/ VN, (20) 

with range measurement a in meters. Note that this result compares favorably with the cruder analytic 
result presented by expression (17). 
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An estimate of how well rj could be estimated from existing ranging data, as discussed at the end of 
Section II, yields: 

a v « 0.00055 - 0.009 (Mariner 9, Viking, Phobos), (21) 

which is slightly smaller than the analytic result given in eq.(18). 

While numerical integration is expected to be more accurate than analytic approximations, one might 
wish to gain more understanding of the planetary SEP effect by doing a realistic analytic calculation 
which improves on the "circular orbit" approximation of Section II. It is natural to eliminate the 
largest sources of error of that approximation. Mars has an orbital eccentricity of .093, and there 
is no fundamental barrier to using elliptical reference orbits. The analytic calculation of Section II was 
redone with elliptical reference orbits for Earth and Mars, working to first order in the eccentricity. 
We also used a more sensible way to include the eccentricity corrections, the method of variation of 
parameters (Robertson &; Noonan 1968), and were able to solve the variation of parameters equations 
for the perturbed orbits of Earth and Mars to fourth order in the eccentricity. Figure 2 shows plots 
of the partial derivative of the Earth-Mars range with respect to the SEP acceleration. Comparison 
of the three curves shows that the eccentricity correction plays a more fairly significant role, than 
one might expect. One reason for this is that the eccentricity corrections turn out to include more 
"secular" matrix elements which are proportional to the time t. Such elements dominate at large times, 
and the eccentricity corrections thereby qualitatively change the nature of the solution in the linear 
approximation. 

— * 

4 Numerical Integration With the Tidal Term A ti d. 

Jupiter's mass jij needs treatment as an adjastable parameter to be fit with the ranging data. This 
is because the octopolar tide of Jupiter acting on the orbits of Earth and Mars produces polarizations 
similar to those produced by the SEP effect, but fortunately having a different Earth-Mars ratio and 
therefore separable from the desired SEP effect. If Jupiter's mass were uncertain by 4 parts in 10 8 , 
its tidal polarization of Mars' orbit would be uncertain by that would be produced by an rj ~ .001, 
for example. But Jupiter's mass is only known to a part in a million, so we must include ji j as a free 
parameter with its own partial. 

In this Section we outline the most accurate calculation. The full equations (8) (and analogue for 
Mars) were numerically integrated including the tidal term And, previously was neglected in Sections 
II and III. To the parameters r Eo ,r Mo ,PEr ,PMr ,PEe ,PMe , &E ~ Qm ,Hs and g v , then, we add jij in 
the covariance analysis, otherwise, the analysis is identical to that of Section III. 

The results are shown in Figure 1, with the solid curve giving the result for from N ranging 
measurements, each with error a meters. For a mission time of the order ten years we find 

a v ~ .0039^/ ViV. (22) 

The estimation of how well rj can be determined from existing ranging data (as discussed in Section 
II) is: 

a v « 0.0012 - 0.02 (Mariner 9, Viking, Phobos), (23) 

about double the result from Section III. 

The covariance analysis gives the expected formal error in jij as well, analogous to eq.(16), with the 
result shown in Fig. 3. For a mission time of order ten years we find a^j ~ 5.7a/-\/N in km^s~ 2 , 
where N, as before, is the number of daily ranging measurements taken during the mission. For 
a = 7.9m, a^j falls below the present accuracy determined from Pioneer 10,11 and Voyager 1,2, namely, 
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Lij = 100fcm 3 /s 2 (Campbell & Synnott 1985) - within two years. This surprising result suggests that 
Earth-Mars ranging can provide an extremely accurate value for the mass of Jupiter, independent of 
the any determination of a SEP effect. 

Following Section II, in order to obtain the realistic estimates for a v and a^j one multiplyes these 
results by y/N /3, giving values of the upper bound of realistic errors for r\ and iij as a v ~ .0012(7 and 
a^j ~ 1.9a km s s~ 2 correspondingly. With mission duration of order ten years this result gives the 
interval for the uncertainties a v and u MJ during the future Mars missions: 

a v « (0.0001 - 0.0012) a, 

» (0.09 - 1.9) a W/s 2 . (24) 

5 Discussion 

The planet Mars has become an object of intensive investigation by many scientists around the world. 
The next flight opportunity during 1996-97 will mark the initiation of a number of new space missions 
to that planet from which we expect to obtain a rich set of data, including spacecraft tracking and 
planetary radar measurements, and allowing precise relativistic celestial mechanics experiments. 

Anticipating these events, we have analyzed the ability for testing SEP violation with Earth-Mars 
ranging. The expected accuracy of the future ranging experiments would put significant constraints on 
theoretical models, including a possible inequality of the Sun's inertial and gravitational masses. Using 
analytic and numerical methods we have shown from covariance analysis that Earth-Mars ranging 
can provide a quality estimate of 77. Indeed, for N ranging measurements with an accuracy of a 
meters, the SEP parameter rj according to covariance analysis can be determined within the accuracy 
a v ~ 3.9 x 10~ 3 a/s/N. The "realistic" estimate for a v based on a "modified worse case" analysis sets 
a conservative limit on the accuracy and indicates that even in unfavorable cases the Sun-Earth-Mars- 
Jupiter system allows for a sensitive test of the Strong Equivalence Principle, qualitatively different 
from that provided by LLR. The mass of Jupiter, fij, can be determined more accurately from a few 
years of Earth-Mars ranging than from Pioneer 10,11 and Voyager 1,2 combined. This analysis shows a 
rich opportunity for obtaining new scientific results from the the program of ranging measurements to 
Mars. 

Efforts are underway at JPL to determine rj from the Mariner 9, Viking and Phobos ranging data. 
This research will modify the theoretical model to include effects due to Saturn. We will perform 
the numerical experiments with combined data collected from the planetary missions, LLR and VLBI. 
And, as the data do not presently include any direct ranging to the Sun, it will be interesting to 
include ranging results to the spacecrafts of the joint US-Russian Solar probe missions scheduled for 
launch in the year 2001, or shortly thereafter. A preliminary determination from combined solar-system 
data, including Mars ranging and lunar-laser ranging, has been reported at a Division of Dynamical 
Astronomy meeting by Chandler et al. (1994). However, we have found that the inclusion of the 
acceleration of eq.(7) in the JPL planetary ephemerides system as postulated (Standish et al. 1993), 
has not been straightforward. The total SEP range signal in Earth-Mars ranging is so complex and 
unique, one should be very conservative about physical interpretation of the results obtained. The 
full scale research of this important experiment is currently underway at JPL. The results reported 
here provide insights into what is being measured, and hence they minimize the possibility of error in 
implementing the planetary SEP effect in the complicated software system. And by concentrating on 
covariance analysis, information needed for the planning of gravitational experiments on future Mars 
missions is obtained. We intend that this paper serve as one of a collection of quides for scientific goals 
and priorities on Mars missions over the next decades. 
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Finally we mention that the analysis of solar ranging data might provide the opportunity for an- 
other fundamental test, namely a Solar system search for dark matter (Nordtvedt 1994; Braginsky 
1994; Nordtvedt 1995). Suppose that dark matter weakly interacts with ordinary matter in a manner 
depending on the specific properties of the matter composing the bodies. The Sun, having an internal 
structure and matter composition which is considerably different from the rest of the inner bodies in 
the Solar system, might then have a different coupling to dark matter, and a corresponding anomalous 
cosmic force F c = msA c acting on the Sun, would produce extra terms in the heliocentric equations 
of motion for the planets - like A^, but fixed in direction, in equation (8). Interesting limits on the 
size of any SEP violation for extended bodies in the Solar system falling toward dark matter may be 
obtainable. This research will be the subject of a subsequent publication. 

We are indebted to our colleagues John Armstrong, Eunice Lau and Skip Newhall for many useful and 
stimulating conversations. MG acknowledges the partial support by an AWU-JPL sabbatical fellowship. 
KLN was supported in part by National Aeronautic and Space Administration throughout Contract 
NASW-4840. SGT was supported by National Research Council, Resident Research Associateship 
award. This work was carried out in part at the Jet Propulsion Laboratory, California Institute of 
Technology, under a contract with National Aeronautic and Space Administration. 
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7 Figure Legends 

Fig. 1. Variation of the uncertainty in the SEP parameter r\ with N, assuming N Earth-Mars range 
measurements, each with uncertainty a meters. The thin dashed curve is the result of the analytic 
approximations of Section II. The thick dashed curve is the numerical integration result of Section III. 
The solid curve comes from the numerical integration described in Section IV where the term A tic i was 
restored and [ij was included in the covariance analysis. 

Fig. 2. Investigation of first order eccentricity corrections (thick dashed curve) to the linear ap- 
proximation of Section II (thin dashed curve). Both of these curves should be compared to the more 
accurate numerical results of Section III (solid curve). All three of these calculations used eq.(8) without 
the And term. The initial conditions taken were different from those in figures 1 and 3. 

Fig. 3. Plot of the uncertainty in the mass of Jupiter versus N from the numerical integration of 
equation (8) plus covariance analysis. See Section IV. 
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